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Abstract: 



A measurement of the proton structure function F^ix, Q"^) is reported for mo- 
mentum transfers squared Q"^ between 4.5 GeV^ and 1600 GeV^ and for Bjorken 
X between 1.8 • 10"'* and 0.13 using data collected by the HERA experiment HI in 
1993. It is observed that increases significantly with decreasing x, confirming 
our previous measurement made with one tenth of the data available in this analy- 
sis. The Q"^ dependence is approximately logarithmic over the full kinematic range 
covered. The subsample of deep inelastic events with a large pseudo-rapidity gap 
in the hadronic energy fiow close to the proton remnant is used to measure the 
"diffractive" contribution to F^. 
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1 Introduction 



The measurement of the inclusive deep inelastic lepton-proton scattering cross section has 
been of great importance for the understanding of quark-gluon substructure of the proton [1]. 
Experiments at HERA extend the previously accessible kinematic range up to very large 
squared momentum transfers, Q"^ > 5-10'* GeV^, and to very small values of Bjorken x < 10"'*. 
In 1993 the first observation was reported of a rise of the proton structure function F^ix, Q"^) 
at low X < 10~^ with decreasing x [2, 3]. Such a behaviour is qualitatively expected in 
the double leading log limit of Quantum Chromodynamics [4]. It is, however, not clarified 
whether the linear QCD evolution equations, as the conventional DGLAP evolution [5] in 
InQ^ and/or the BFKL evolution [6] in ln(l/a;), describe the rise of or whether there is 
a significant effect due to nonlinear parton recombination [7]. The quantitative investigation 
of the quark-gluon interaction dynamics at low x is one of the major challenges at HERA. 
It requires high precision for the F^ measurement and an investigation of the hadronic final 
state. 

The structure function F^ix^Q"^) is derived from the inclusive lepton-proton scattering 
cross section. It depends on the squared four momentum transfer Q"^ and the scaling variable 
X. These variables are related to the inelasticity parameter y and to the total squared centre 
of mass energy of the collision s as = xys with s = iEgEp. In 1993 the incident electron 
energy was Eg = 26.7 GeV and the proton energy was Ep = 820 GeV. A salient feature of 
the HERA collider experiments is the possibility of measuring not only the scattered electron 
but also the complete hadronic final state, apart from losses near the beam pipe. This means 
that the kinematical variables x, y and Q'^ can be determined with complementary methods 
which experimentally are sensitive to different systematic effects. The comparison of the 
results obtained with different methods improves the accuracy of the F2 measurement. A 
convenient combination of the results ensures maximum coverage of the available kinematic 
range. 

In this paper an analysis is presented of inclusive deep inelastic scattering (DIS) data 
taken by the HI collaboration in 1993 with an integrated luminosity of 0.271 pb~^. During 
1993, its second year of operation, HERA delivered an integrated luminosity of an order of 
magnitude larger than in 1992. With these increased statistics the accessible kinematic range 
has been extended and the behaviour of F2 has been investigated at a new level of precision. 
A similar measurement was published recently by the ZEUS collaboration [8]. 

The structure function results presented here are more precise with a typical systematic 
error of 10% than the previous HI data. The following new information is provided: i) 
analyzing data with shifted vertex position along the proton beam line, the first DIS data for 
Q'^ between 5 and 10 GeV^ at a; < 0.001 is obtained in a region where new effects could be 
observed; ii) the analysis extends down to y ~ 0.01 which allows for the first time to approach 
with HERA data the region of the fixed target lepton-proton scattering experiments; iii) due 
to the increased statistics, the first precise HI measurement of F2 beyond Q"^ ~ 100 GeV^ 
can be presented. An analysis is presented of the dependence of F2 on the effective mass W 
of the virtual photon-proton system which, due to the very large energy s at HERA, can be 
performed in a new kinematic range. 

A first measurement is given of the diffractive contribution to F2 by analyzing the sub- 
sample of about 6% of the DIS events which exhibit no activity in the forward detector region 
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where "forward" denotes the direction of the proton beam. Characteristics of the events with 
a large pseudo-rapidity gap in the hadron flow close to the proton remnant direction have 
already been studied at HERA [9, fO]. 

The paper is organized as follows. After a brief introduction to the kinematics of the 
inclusive ep scattering process (sec. 2), the HI apparatus and the 1993 data taking (sec. 3) 
are described. Then the event selection including the background rejection (sec. 4) and the 
detector response (sec. 5) are discussed. Section 6 describes the F2 analyses and discusses the 
results. In section 7 the diffractive data analysis is presented. The paper is summarized in 
section 8. 



2 Kinematics and Analysis Methods 

The kinematic variables of the inclusive scattering process ep eX can be reconstructed in 
different ways using measured quantities from the hadronic final state and from the scattered 
electron. The choice of the reconstruction method for Q'^ and y determines the size of 
systematic errors, acceptance and radiative corrections. The basic formulae for Q"^ and y 
used in the different methods are summarized below, x being obtained from Q'^ = xys. For 
the electron ("E") method 

. . 2^e ^2 ^I'sin^^, 



where the electron polar angle is defined with respect to the incident proton beam direction 
(+z axis). The resolution in Ql is 4% while the precision of the yg measurement degrades 
as 1/ye- Thus the electron method cannot be used for yg < 0.05. In the low y region it 
is, however, possible to use the hadronic methods for which it is convenient to define the 
following variables 

S = ^ (E, - p,,) (p!^r = (J2p.,ir + C£Py,if- (2) 

i i i 

Here E,Px,Py,P;, are the four-momentum vector components of each particle and the sum- 
mations extend over all hadronic final state particles. The standard definitions for y^ [11] 
and 9h are 

S Oh T. 

Vh = tan — = —r. (3) 

^ 2E, 2 ^ ' 

The combination oiy^ and Ql defines the mixed method [12] which is well suited for medium 
and low y measurements. It was used in the previous HI analysis of F^ [2]. The double-angle 
("DA") method [13] makes use only of 9g and with 



tan(g,/2) ^ coijeji) 

tan(^j2) + tan(^ft/2) nan(^j2) + tan(^ft/2) 



The method is rather insensitive to the absolute energy calibration of the detector. It has 
good resolution at large Q"^ where the jet energies are high but the resolution deteriorates 
for X < 0.001. The formulae for the S method [14] are constructed requiring Q'^ and y 
to be independent of the incident electron energy. Replacing in eq.3 the quantity 2Ee by 
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S + E'g(l — cosOg), as allowed by the conservation of the total E — of the event^, and 
(1 — Ue) by (1 — ys) in eq.l, one obtains: 

S + - cos6'e) 1 - 

For non radiative events and are equivalent at low y. However, the modified quantity 
is less sensitive than to the imprecise hadronic measurement at high y since the S term 
dominates the E — of the event. Therefore the S method can be applied over the full 
kinematic range considered in this paper. 

All these methods were utilized to measure F^. Resolution values of the x and Q"^ variables 
reconstructed with the HI detector are discussed for all methods at the end of section 5. 



3 The HI Detector and the Data Taking in 1993 
3.1 The HI Detector 

The HI detector [15] is a nearly hermetic multi-purpose apparatus built to investigate the 
inelastic high-energy interactions of electrons and protons at HERA. The structure function 
measurement relies essentially on the inner tracking chamber system, on the backward elec- 
tromagnetic calorimeter, and on the liquid argon calorimeters which are described briefiy 
below. For the luminosity measurement see section 3.2. 

The tracking system includes the central tracking chambers, the forward tracker modules 
and a backward proportional chamber. These chambers are placed around the beam pipe at 
z positions between —1.5 m and 2.5 m. A superconducting solenoid surrounding both the 
tracking system and the liquid argon calorimeter provides a uniform magnetic field of 1.15 T. 

The central jet chamber (CJC) consists of two concentric drift chambers covering a polar 
angle of 15° to 165°. Tracks crossing the CJC are measured with a transverse momentum 
resolution of Spx/pr < 0.01 ■ pT/GeV. The CJC is supplemented by two cylindrical drift 
chambers at radii of 18 and 47 cm, respectively, to determine the z coordinate of the tracks. 
To each of the z drift chambers a proportional chamber is attached for triggering. The inner 
one (CIP) was used in addition to estimate residual photoproduction background. 

A tracking chamber system made of three identical modules measures charged particles 
emitted in forward direction (3° to 20°). This forward tracker (FT) is used to determine the 
vertex for the events which leave no track in the CJC. This extends the analysis into the 
larger x region. 

In the backward region a four plane multiwire proportional chamber (BPC) provides a 
space point for charged particles with a spatial resolution of about 1.5 mm in the transverse 
plane. The polar angle acceptance of the BPC ranges from 155° to 174.5°. The reconstructed 
space point together with the z vertex position defines the polar angle of the scattered elec- 
trons. 

The BPC is attached to the backward electromagnetic calorimeter (BEMC) which is made 
of 88 lead/scintillator stacks with a size of 16 X 16 cm^ and a depth of 22 radiation lengths, 

"'defined as (E — Pz) = ^ (Ej — Pz,j), the sum extending over all particles j of the event. 
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corresponding to about one hadronic interaction length. The angular coverage of the BEMC 
is 155° < 9 < 176°. A 1.5 cm spatial resolution of the lateral shower position is achieved 
using four photodiodes which detect the wavelength shifted light from each of the scintillator 
stacks. A scintillator hodoscope (TOF) situated behind the BEMC is used to veto proton- 
induced background events based on their early time of arrival compared with nominal ep 
collisions. 

Hadronic final state energies and electrons at high Q'^ are measured in the highly seg- 
mented liquid argon (LAr) calorimeter [16] which covers an angular region between 3° and 
155°. The LAr calorimeter consists of an electromagnetic (e.m.) section with lead absorber 
and a hadronic section with stainless steel absorber. The electromagnetic part has a depth 
between 20 and 30 radiation lengths. The total depth of both calorimeters varies between 
4.5 and 8 interaction lengths. The most backward part of the LAr calorimeter is a smaller 
electromagnetic calorimeter (BBE) ^ which covers the polar angle range from 146° to 155°. 

The DIS events in the backward region (low Q"^ data, Og > 155°) were triggered by an 
energy cluster in the BEMC (E > 4 GeV) which was not vetoed by an out of time signal in the 
TOF. The high Q"^ events were triggered by requiring an e.m. energy cluster with E > 8 GeV 
in the LAr calorimeter. For lower energy thresholds (> 5 GeV) the events were also triggered 
if there was simultaneously a tracking trigger. The trigger efficiency has been determined 
from the data using the redundancy of the HI trigger system which relies on calorimetry and 
tracking. In the region of the final F2 data presented below, i.e. for E'^^ > 10.6 GeV, the DIS 
electron trigger efficiency is 100% within the errors. 

3.2 Data Samples and Luminosity 

In 1993 HERA was operated with 84 colliding electron and proton bunches. Sixteen pilot 
bunches, 6 proton and 10 electron, underwent no ep collision. From these the beam induced 
background could be estimated. A small part of the data was taken with the nominal inter- 
action position shifted in z by -|-80 cm in order to reach Q"^ values as low as 4 GeV^. The 
lower Q"^ region was covered also by analyzing events which originated from the so-called 
early "proton satellite" bunch which collided with an electron bunch at 2; ~ +62 cm. The 
satellite bunch data amount to ~ 3% of the total data. Both event samples with shifted z 
vertex positions were analyzed independently and the results combined. 

To reduce the systematic errors of the F2 measurement, a strict data selection was per- 
formed based on the behaviour of the main detector components. In particular, data taken 
during a period in which there was no magnetic field due to a failure of the coil were excluded 
(~ 0.15 pb~^). The number of accepted events per luminosity was checked to be constant 
within statistical errors during the full data taking period. 

The luminosity was determined from the measured cross section of the Bethe Heitler 
reaction e~p e~prj. The final state electron and the photon are simultaneously detected 

^The HI detector calorimetry and electron detection is thus split into two parts. At large angles, d > 155°, 
the electron is measured in the BEMC and preceding tracking chambers. This defines for the subsequent 
analysis the low data sample. Electrons at lower polar angles, 6 < 155°, are detected in the LAr calorimeter 
and the forward and central chambers defining the high event sample. For the structure function analysis 
both data samples were combined but the analyses of the low and high data required partially specific 
techniques and cuts. Special attention had to be paid to the BEMC-BBE transition region where the electron 
energy is often shared between the two calorimeters. 
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in small electromagnetic calorimeters (electron and photon "taggers") close to the beam 
pipe but at a large distance from the main detector (at z = —33 m and z = —102 m). 
This coincidence method gives the best on-line luminosity estimate but it is sensitive to 
variations in the electron beam optics through the electron tagger acceptance. Therefore the 
final luminosity determination was performed based on the hard photon bremsstrahlung data 
using the photon tagger only. The results of these two methods differ by 1.1%. This is well 
below the overall systematic error which is estimated to be equal to 4.5% [17]. 

An independent check of the luminosity measurement has been performed by determining 
the cross section for QED Compton events which leave a scattered electron and a photon in 
the backward part of the apparatus. For the 1993 data this yields a luminosity value which 
is 4% lower than the one from the Bethe Heitler cross-section measurement, with a combined 
systematic and statistical error of 8% [17, 18]. 

The integrated luminosity for the nominal vertex data used in this analysis is 0.271 pb~^, 
and 2.5 nb~^ for the special data taking with the shifted vertex position. The luminosity of 
the satellite data was obtained from the measured luminosity for the nominal vertex data 
multiplied by the efficiency corrected event ratio in a kinematic region common to both data 
sets. The error of that luminosity determination was estimated to be 8.5%. 

4 Event Selection and Monte Carlo Simulation 

4.1 Event Selection 

The event selection criteria can be divided into four categories: i) electron identification, 
ii) event vertex requirement, iii) kinematical constraints and iv) background rejection. The 
latter is discussed in the subsequent section. A summary of the selection cuts is given in 
tables 1 and 2. 

i) For the electron identification in the backward region of the HI detector, for Og > 
155°, the most energetic cluster in the BEMC is selected. Its center of gravity is required 
to be at a radial distance ("CLBP") smaller than 4 cm from a reconstructed BPC point. 
The lateral size ("ECRA") of that cluster, calculated as the energy weighted radial distance 
of the cells from the cluster centre, has to be smaller than 4 cm. In the LAr calorimeter, 
for 9e < 155°, several electron identification algorithms have been developed which have 
the common feature of demanding significant electromagnetic energy deposited in a confined 
region (table 2). The E method required an isolated electromagnetic cluster with 64 > 50% 
and a minimum relative energy ("£3") of 3% in the first e.m. layer of the LAr calorimeter 
which amounts to about three radiation lengths [19]. Here 64 is the energy sum over the four 
most energetic cells of a cluster divided by its total energy. For the S analysis an e.m. cluster 
with 64 > 65% was considered to be an electron. 

In the transition region between the BEMC and the BBE, for Og about 155° ± 2°, a 
dedicated E and DA analysis was performed in which the electron finding was based on 
topological criteria only, i.e. without making use of E'^^ explicitely. For the S analysis a 
fiducial cut was applied if the electron cluster was shared between the BEMC and the BBE. 
Because of the smearing of the z vertex position there is no loss in the x, Q"^ coverage, only 
a reduction in statistics. 
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The electron identification efficiency for DIS events, as determined from Monte Carlo 
simulation, is better than 97% apart from the BEMC-BBE transition region where it amounts 
to approximately 92%. The scattered electron is considered to be the cluster with the largest 
energy. If two electron candidates are present in the same event, the one of highest energy 
is selected. This introduces a misidentification probability of at most 3% [20] at the smallest 
electron energy which was considered in the systematic error calculation. 



ii) The event vertex position is needed for a precise determination of the kinematics. It 
was defined by at least one well measured track in the CJC or in the FT crossing the beam 
axis. The z vertex region for the events with interactions in the nominal time intervals was 
z = ( — 5 ± 30) cm. For the satellite bunch data the requirement was z = (+62 ± 20) cm 
and for the shifted vertex data z = ( + 75 ± 25) cm. The vertex reconstruction efficiency, as 
determined from the data, is 97% at y > 0.2. It decreases smoothly to an average of 73% 
in the lowest y bins (y ~ 0.01) since with decreasing y the hadronic particles are emitted at 
smaller polar angles. The vertex reconstruction efficiency for the shifted vertex data is only 
about 10% lower than for the nominal vertex data since 2; ~ 75 cm is still inside the central 
tracking system. 



iii) The basic kinematic constraints were a maximum electron scattering angle of 173° 
(175.2° for the shifted vertex data) and minimum energy requirements in order to ensure 
high trigger efficiency and a small photoproduction background. For the high Q"^ data, the 
minimum energy requirement was replaced by introducing a cut rejecting events at high yg. 
At low Q"^, the minimum electron candidate energy E'^^ was 10.6 GeV for the E analysis and 
8 GeV for the S analysis. That difference is an example of selection criteria differences arising 
due to the combination of several complete and independent analyses. The S analysis also 
required a total missing transverse momentum smaller than 15 GeV and a total E — P^ between 
30 and 75 GeV. In analyzing the BEMC-BBE transition region alternative requirements, like 
^ogiyoAl Vh) < 0.5, were introduced avoiding the use of E'^^ in the selection. The latter 
conditions rejected photoproduction background events and DIS events with an energetic 
photon radiated along the electron beam direction, thus reducing the radiative corrections to 
F2. 





E method 


S method 


DA method 




< 173 


< 173 


< 173 


E'JGeV 


> 10.6 


> 8 


> 10.6 


ertex /cill 


-5+ 30 


-5+ 30 


-5+ 30 


ECRA/cm 


< 4 


< 4 


< 4 


CLBP/cm 


< 4 


< 4 


< 4 


(^-P,)/GeV 




< 75, > 30 




Oh/" 






> 40 


selected events 


22500 


24100 


22200 



Table 1: Summary of event selection criteria for the low Q"^ data, i.e. for 9g > 155°. For the 
dehnitions of abbreviations see the text. 
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E method 


S method 


DA method 




< 150 


< 155 


< 150 


Ve 


< 0.6 


< 0.8 


< 0.8 


ertex /cill 


-5± 30 


-5± 30 


> 1 track 


electron ident. 


64 > 50% 


64 > 65% 


min Rci 




63 > 3% 


> 8 GeV 


Ei>8 GeV 


Pt/Pt 


> 0.2 


— 


— 


iE-P;)/Ge\ 


— 


> 30 


> 30 


missing px/GeV 




< 15 


< 10 - 25 


selected events 


880 


1100 


1038 



Table 2: Summary of event selection criteria for the high Q"^ data. Further topological 
requirements based on the LAr calorimeter were imposed to reject background. The quantity 
Rci in the DA method denotes the average over calorimeter cell energies of the polar angle 
distance between a cell and the cluster centre which is a measure of the lateral extension of 
the calorimeter energy deposition. For further definitions see the text. 

4.2 Background Rejection 

4.2.1 Non-ep Background 

At low Q"^ the main sources of non-ep background are due to proton beam interactions with 
residual gas and beam line elements upstream of the HI detector. Beam wall events in the 
interaction region were rejected by requiring the reconstructed vertex to be centered on the 
beam axis. An efficient reduction of the remaining background is provided by the minimum 
energy and the vertex requirements discussed above. In addition two algorithms have been 
developed based on the central tracking information. An event was rejected either if the 
fraction of tracks which did not point to the reconstructed vertex was too high, or if at least 
two tracks pointed to the backward direction and crossed the beam line at z < —50 cm. Both 
algorithms gave similar results. They introduced a loss of a few per cent of genuine DIS events 
which was controlled by eye scanning and corrected for. From the study of pilot bunches, the 
residual background was estimated to represent less than 1% of the total number of selected 
events. 

At high Q'^ the main background is due to muons travelling off axis parallel to the proton 
beam. These are produced by proton beam halo interactions and occasionally generate an 
electromagnetic shower in the LAr calorimeter. Requiring a reconstructed vertex rejects most 
of them. Further rejection was obtained by requiring more than 10 GeV deposited outside the 
calorimeter region which contains the electron candidate cluster. Residual cosmic ray event 
candidates were rejected with similar topological requirements. The remaining background 
was estimated to be smaller than 2% by a visual event scan and also by analyzing the pilot 
bunch data. 

4.2.2 Photoproduction Background 

The only significant background to DIS from ep interactions is due to photoproduction events 
at ~ where the scattered electron escapes the detector along the beam pipe but in which 
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an energy cluster from the hadronic final state fakes an electron. About 10% of these events 
are identified as photoproduction background if the scattered electron is found in the electron 
tagger. 




Figure 1: Distributions for tagged photoproduction events: a) Energy distribution of the 
scattered electron energy in the electron tagger, b) and c) energy and angle distributions of 
the fake electron candidate in the BEMC. The statistical errors of the Monte Carlo simulation 
are of the same order of magnitude as the experimental ones. 

To estimate this background, photoproduction events were simulated corresponding to 
the luminosity of the data. The "soft" vector meson contribution was simulated using the 
RAYVDM [21] program, and the "hard" scattering part using the PYTHIA [22] program. 
For hard scattering interactions, direct and resolved processes and the production of heavy 
quarks were included. The relative contributions of both were adjusted to agree with the total 
photoproduction cross section analysis [23]. The simulation both in shape and normaliza- 
tion [24, 25] is in good agreement with the tagged photoproduction data. This is illustrated 
in fig.l showing the scattered electron energy distribution in the tagger and the energy and 
angle distributions of the fake electron detected in the BEMC. The photoproduction back- 
ground was subtracted statistically bin by bin. The highest contamination is in the lowest 
Q"^, X bin and amounts to (9 ± 4)% for E'^ > 10.6 GeV. Only three bins have a contamination 
larger than 3% [26]. 

An independent analysis [27] made use of the tracking information to identify that part 
of the photoproduction background which originates from photons in a low Q"^ interaction. 
These can mimic a DIS electron due to conversion in the CJC end fiange. Such events were 
rejected by requiring a hit in the innermost tracking chamber which is crossed by the particle 
before reaching the CJC. The remaining background from charged pions and conversions in 
the beampipe was estimated based on a comparison of the tagged photoproduction events 
and the Monte Carlo simulation which allowed to subtract this background part statistically. 
This background estimation gave consistent results with the Monte Carlo method described 
above. 
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4.3 Monte Carlo Simulation 



More than half a million Monte Carlo events were generated using D JANGO [28] and different 
quark distribution parametrizations, corresponding to an integrated luminosity of approxi- 
mately 1.5 pb~^. The program is based on HERACLES [29] for the electroweak interaction 
and on LEPTO [30] to simulate the hadronic final state. HERACLES includes first or- 
der radiative corrections, the simulation of real bremsstrahlung photons and the longitudinal 
structure function. The acceptance corrections were performed using the MRSH parametriza- 
tion [31], which is constrained to the HERA F2 results of 1992. LEPTO uses the colour dipole 
model (CDM) as implemented in ARIADNE [32] which is in good agreement with data on the 
energy fiow and other characteristics of the final state as measured by HI [33] and ZEUS [34]. 
For the estimation of systematic errors connected with the topology of the hadronic final 
state, the HERWIG model [35] was used in a dedicated analysis. Based on the GEANT 
program [36] the detector response was simulated in detail. After this step the Monte Carlo 
events were subject to the same reconstruction and analysis chain as the real data. 

5 Calibration and Kinematical Distributions 

The structure function measurement requires an understanding of the reconstruction of the 
energies and angles of the scattered electron and of the hadronic final state because all are 
used to define the kinematics of the interaction. Figs. 2 and 3 display the distributions 
of the reconstructed energy E'^^ and angle Og of the scattered electron and of the hadronic 
variables and 9^. Fig. 2 shows the distributions for the high statistics, low Q"^ event sample 
with the electron measured in the BEMC. The small contamination due to photoproduction 
background is also displayed. The Monte Carlo distributions are normalized to the luminosity 
measurement. Agreement between data and Monte Carlo simulation was obtained after 
an iteration of the structure function parametrization (see below). The remaining small 
discrepancies, as visible for example in the distribution of around 0.3, have negligible 
infiuence on the acceptance calculation for F2. Fig. 3 shows the same distributions for the 
smaller sample at large Q'^ with the electron detected in the LAr calorimeter. Data and 
simulation agree well. 

A crucial part of the F2 analysis is the absolute energy calibration. The determination 
of the energy scale of the BEMC is based on the observed shape of the kinematic peak, see 
fig. 2a, for energies between 22 and 28 GeV. Taking into account stack to stack variations of 
the response and correcting for dead material in front of the BEMC as well as for energy losses 
between the stacks, the absolute scale of the E'^^ measurement in the BEMC was determined 
with 1.7% systematic and 1% statistical accuracy [37]. The result was cross checked with 
the double angle method to reconstruct the electron energy and agreement was found to 
within the quoted uncertainty. For the high Q"^ data the double angle method was used to 
refine the energy scale determined by test beam measurements [38]. The resulting systematic 
uncertainties in E'^^ are smaller than 5% for Og between 150° and 155° (BBE region) and 3% 
for electrons in the barrel part of the calorimeter (Og < 150°). 

The hadronic energy scale in the liquid argon calorimeter is presently known to 6% as 
determined from studies of the transverse momentum balance of DIS events. Test beam 
data of pions between 3.7 GeV and 205 GeV showed agreement on the 3% level with the 



12 





Figure 2: Distributions of the lunematical quantities E'g,9e,y-E and 9^ for the low Q"^ data 
(closed circles) and the Monte Carlo simulation (open histogram). The Monte Carlo cal- 
culation is normalized to the luminosity. The hatched histogram is the estimation of the 
background due to photoproduction processes. 
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Figure 3: Distributions of the lunematical quantities El,9e,y-E and 9^ for the high Q"^ data 
(9g < (closed circles) and the Monte Carlo simulation (open histogram) using the MRSH 

parametrization [31] as input structure function. For the energy distribution Q'j > 250 GeV^ 
is required. For the 9^ distribution a cut at > 0.05 is applied to exclude the region where 



is not measured. 
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Monte Carlo description [39]. The calibration parameters were determined from Monte Carlo 
simulated jet data [15, 40]. The calorimetric energy due to charged particles in the central 
detector region was replaced by the more precise momentum measurement of the C JC which 
renders sensitive to charged hadrons of low momentum [41]. This improved the resolution 
of Vh by about 40% at low Q"^ and high y. 

An extension to higher values from 0.4, as in the previous HI analysis [2], to 0.7 was 
achieved by using the S method. The behaviour of can be studied with the high y data 
by means of the ratio ys/ye (fig. 4a) using the excellent resolution of yg (— 4%) as a reference. 
As can be seen, both the mean value and the measured resolution (13%) are well reproduced 
by the Monte Carlo simulation. Good agreement is also achieved in the p^/Pt distribution 
at high y (fig. 4b) although the hadronic final state particles are of rather low energy and on 
average are emitted into the backward direction. The resolution of (~ 35%) is wider than 
the resolution. The coverage of the full y range necessitates an understanding of the effect 
of the calorimeter noise on F2. For low y this is sensitive in particular in the backward part of 
the detector where any energy deposition of 300 MeV generates an additive contribution to 
y of about 0.01. Fig. 2c suggests that y^ values below 0.01 can be reconstructed. In fig. 4c the 
y's/Vgen distribution is shown for 0.005 < ys < 0.015 including and excluding the calorimeter 
noise contribution for the reconstruction of y. The distributions have a maximum near 1 and 
the resolution obtained is 27% when the noise is included. Note that for low y the quantity 
ye cannot be considered as a measure for the true y since its resolution function is distorted. 
The contribution of the noise which has been obtained from experimental data is apparently 
small. The tail of the distribution at large ys/ygen is related to the finite granularity of the 
calorimeter [42]. Good agreement between data and Monte Carlo in the low y region has 
been achieved. This can be seen, for example, in the distribution of the ratio p!^/p^, fig.4d. 
Therefore, the F2 measurement could be extended to y = 0.01 and, for the first time, the 
X region measured at HERA reaches the domain of larger x values covered in fixed target 
experiments [43, 44]. 

The electron angle Og has been determined for the lower Q'^ data from the vertex position, 
reconstructed with the central and forward chambers, and the hit in the BPC closest to 
the position of the electromagnetic cluster in the BEMC. The precision of the z vertex 
measurement for most of the events is about 1 cm. The Og measurement, in part of the 
kinematic range, could be validated and cross calibrated with the innermost tracking chamber 
and the cluster position [27]. Potential shifts of the Og values are estimated to be smaller than 
2 mrad and a resolution of 2.5 mrad is achieved. For the high Q'^ data, the vertex position and 
the cluster centre define Og within 5 mrad accuracy and 7 mrad resolution as cross checked 
with the central jet chamber. For the intermediate Q"^ analysis, the 9g measurement is made 
with the CJC and the z drift chambers. 

The hadronic angle 9^ is reconstructed according to eq.3 from the energy deposited in 
the calorimeter cells. The 9^ Monte Carlo distribution depends crucially on the shape of 
the input structure function. Agreement between data and Monte Carlo was obtained after 
iteration of the input by starting with the MRSH distribution, extracting and fitting F2, and 
reweighting the Monte Carlo event distributions by the ratio of the fitted F2 to the MRSH 
parametrization. The weights applied deviate from 1 by less than 20%. While this procedure 
significantly improved the description of the event distributions shown in fig. 2, it had less 
than 5% infiuence on the structure function derivation because the dependence of the F2 
analysis on the assumed shape of F2 is due only to smearing corrections which are less than 
30% in the kinematic range and for the binning used here. 
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Figure 4: Experimental and Monte Carlo distributions of a) ys/ye and b) Px/px in the high 
y range for Q| > 10 GeV'^. Monte Carlo distributions of c) ys/ygen with and without (w/o) 
the calorimeter noise included and d) experimental and Monte Carlo distributions ofp^/pfp 
at very low y. The px ratio distributions at low y (d) peak at about 0.8 due to losses of 
hadrons in the beam pipe both in the data and in the Monte Carlo simulation. 
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Detailed studies were performed of the resolution and systematic shifts of the recon- 
structed X and Q'^ values. In table 3 the average resolution values are quoted for the three 







Xe 




Xda 


Qe 






y > 0.15 


low g2 


12 


13 


39 


4 


11 


12 




high 


9 


10 


9 


3 


9 


4 


y < 0.15 


low g2 


40 


19 


23 


4 


7 


3 




high 


37 


16 


12 


3 


5 


4 



Table 3: Relative resolution values in per cent for the three methods used to reconstruct x 
and Q^. 

methods used to determine F2. The binning was adapted to the resolution in x and 

to the statistics in Q"^. A rather fine grid is obtained at low Q"^ with 8 bins per order of 
magnitude for Q'^ between 7.5 and 133 GeV^. Statistics allowed for four more bins, equidis- 
tant in logQ^, between 133 and 2070 GeV^. Independently of Q"^ the binning in x was 6 (4) 
bins per order of magnitude for x < (>)10~^. The F2 measurement is based on kinematic 
variables which in the full (x, Q"^) range are reconstructed with systematic shifts smaller than 
5% and relative resolutions better than 20%, apart from three edge bins at very low y with 
a resolution better than 30%. 



6 Structure Function Measurement 
6.1 Comparison of F2 Analyses 

The structure function F^ix, Q"^) was derived after radiative corrections from the one-photon 
exchange cross section 

{2-2y + -^)F2{x,Q^) = K-F2 (6) 



dxdQ"^ Q^x' 1 + R' 

Effects due to Z boson exchange are smaller than 5% and were corrected for at high Q"^. 
The structure function ratio R = F^/lxFi — 1 has not been measured yet at HERA and was 
calculated using the QCD relation [45] and the MRSH structure function parametrization. 
At lowest X and Q"^ the assumed R values are about 0.3, all values being quoted in the F^ data 
summary table, see below. Compared to the previous HI analysis [2] the F^ measurement 
has been extended to lower and higher Q"^ (from 8.5 — 60 GeV^ to 4.5 — 1600 GeV^) and also 
to larger x. 

The determination of the structure function requires the measured event numbers to be 
converted to the bin averaged cross section based on the Monte Carlo acceptance calcula- 
tion. The mean acceptance ^ was 0.89. All detector efficiencies were determined from the 
data utilizing the redundancy of the apparatus. Apart from very small extra corrections all 
efficiencies are correctly reproduced by the Monte Carlo simulation. The bin averaged cross 



■^The data sample contains a fraction of 6% of large rapidity gap events. With a special rapidity gap 
Monte Carlo program [46] acceptances were calculated and they agree within statistical errors with those 
calculated using the standard DIS simulation. Thus only the DIS Monte Carlo simulation has been used for 
the calculation of the acceptance corrections. 
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section was corrected for higher order QED radiative contributions and a bin size correction 
was performed. This determined the one-photon exchange cross section which according to 
eq.6 led to the values for F^ix^Q"^). 

Several methods were used to derive F^, each with different advantages: i) the E method 
which has the best resolutions on x and at large y and is independent of the hadron 
reconstruction, apart from the vertex requirement; ii) the S method which has small radiative 
corrections and extends from very low to large y values; iii) the DA method which is less 
sensitive to the energy scales. A complete analysis [47] has been performed based on the mixed 
method which agrees very well with the other structure function determinations presented in 
this paper. The application of different methods was important to check that the systematic 
errors were correctly evaluated. 

The F2 analyses for the shifted vertex and the satellite bunch data were performed with 
only the E method. The results were found to be in good agreement with each other and with 
the F2 obtained from the nominal vertex data in the region of overlap. Taking into account 
correlations between the systematic errors, the two data sets were combined. Thus the first 
structure function data for ~ 5 GeV^ and x >2 ■ 10"'* are presented here. 

All structure function values are shown in fig. 5. As can be seen the agreement between 
the analyses is very good. A combination of F.f and Ff was chosen for the final result and 
F2^ displayed to demonstrate consistency. 

6.2 Systematic Errors 

The systematic errors considered are: 

• A potential miscalibration of the electron energy by 1.7% in the BEMC, by 5% in the 
BBE, and by 3% in the central calorimeter region. 

• A 6% scale error for the hadronic energy in the LAr calorimeter, the effect of which is 
reduced due to the joint consideration of tracks and calorimeter cells for the S analysis. 
A 20% scale error was assigned to the energy of the hadronic final state measured in the 
BEMC. These numbers include uncertainties due to the noise treatment for the LAr 
calorimeter and the BEMC. 

• A shift of up to 2 mrad for the electron polar angle in the BEMC region and of at most 
5 mrad in the LAr calorimeter. Errors of the energy and angular resolutions are taken 
into account in addition to possible shifts of mean values. 

• Apart from the electron identification all efficiencies were determined from the data and 
compared with the Monte Carlo simulation. Agreement between the experimental and 
the simulated values for the individual efficiencies (trigger, TOE, vertex, CJC track, 
BPC hit, BEMC and LAr calorimeter cluster reconstruction) was found to be better 
than 2%. An overall error of 4% was assigned due to the imperfect description of the 
various efficiencies. A larger error of at most 12% was added to account for the vertex 
reconstruction efficiency variation at large x. 
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Figure 5: Measurement of the structure function F^ix^Q"^) with the electron (closed circles), 
the S (open circles) and the double angle method (open squares). The inner error bar is 
the statistical error. The full error represents the statistical and systematic errors added in 
quadrature, not taking into account the 4.5 % systematic error of the luminosity measurement. 
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Figure 6: Measurement of the proton structure function F^ix^Q"^). The inner error bar is 
the statistical error. The full error represents the statistical and systematic errors added in 
quadrature not taking into account the 4.5% systematic error on the luminosity measure- 
ment. Open circles and triangles represent NMC and BCDMS measurements, respectively. 
A smooth transition becomes apparent from the NMC and BCDMS data (open circles and 
triangles, respectively) to the HI data. The curves represent a phenomenological ht to all 
data, see text. 
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• An error of up to 3% in the radiative correction due to uncertainties from the hadronic 
corrections, the cross section extrapolation towards = 0, second order corrections 
and the absence of the soft photon exponentiation in the HERACLES Monte Carlo. 
The accuracy was cross checked by comparing the Monte Carlo estimate [48] with 
TERAD [49] and also with TERAD supplemented by a leading log higher order cal- 
culation [50]. A direct estimate has been made comparing in their overlap region the 
cross sections derived from the E and the S methods. Note that most of the Compton 
events do not contribute here due to the vertex requirement. 

• The structure function dependence of the acceptance and bin size corrections which was 
controlled to better than 3%. The comparison of the different simulation models of the 
hadronic final state mentioned above was used to assign an additional 3% systematic 
error to the hadronic methods. 

• The uncertainty due to photoproduction background was assumed to be smaller than 
half the correction applied, i.e. smaller than 5%. This affects only the highest y bins 
at lower Q"^ . 

• Statistical errors in the Monte Carlo acceptance and efficiency calculations were com- 
puted and added quadratically to the systematic error. 

• The shifted vertex data sample required to assign a 15% error for the vertex efficiency 
correction due to lack of experimental statistics. For the higher statistics satellite data 
the additional luminosity uncertainty implies that both data samples with shifted z 
vertex position contribute with about the same systematic errors to the final low Q"^ 
structure function measurement. 

6.3 The Structure Function F2(a;,(52) 

The combination of different methods meant that the largest kinematic range could be covered 
in an optimum way for the measurement of F^ix, Q"^). Fig. 6 presents the combined structure 
function data. Data at lower x in the full Q"^ range are obtained with the E method. For 
y < 0.15 the S method was used. Double counting of the events was reduced to the level of a 
few per cent by introducing to the S analysis in each Q'^ bin an Xg cut near to y = Q"^ / sx ~ 
0.15. The results are in good agreement with the previous HI publication [2]. 

The structure function values are given in table 4 with their statistical and systematic 
errors. For low Q"^ values, apart from the lower statistics shifted vertex data, the systematic 
error of is about 10% and twice as large as the statistical error, while for the larger Q"^ 
data the statistical error dominates. 

Fig. 6 shows that the structure function F^ rises steeply with x decreasing to a; = 1.8 •10"'*. 
As can be seen in fig. 7 the dependence of F^ on Q"^ at fixed x is weak. The violation of 
scaling appears to be stronger for smaller values of a;, a trend already observed in fixed 
target experiments for x < 0.1. Compared to the recent data from the ZEUS experiment, an 
extension of the kinematic range is achieved towards low x and low Q^, see fig. 7. 

The X and Q'^ behaviour of F2 can be described by a phenomenological ansatz of the type 
F,(x, Q'') = [a.x' + C-X'' -(l + e-y^)- (InQ' + /In' Q')] ■ (1 - xy . (7) 
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Table 4: Proton structure function F^ix, Q"^) with statistical and systematic errors. All points 
have an additional scale uncertainty of 4.5% due to the luminosity determination. Q"^ is given 
in GeV^. The values of R result from a calculation based on the QCD prescription using the 
MRSH parton distribution parametrization. 
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10^ 10^ 10^ 10^ 

Figure 7: Measurement of the proton structure function F^ix^Q"^). The HI data (closed 
points) is consistent with the ZEUS result (open squares) but extends further towards low 
X and also to lower Q"^ in different x bins. With a small correction the ZEUS data was 
shifted to the Ell x values by using the parametrization given in [8]. The curves represent a 
phenomenological ht to the HI, NMC and BCDMS data, see text. The F^ values are plotted 
with all but normalization errors in a linear scale adding a term c{x) = Q.Q{ix — 0.4) to F^ 
where is the bin number starting at = 1 for x = 0.13. 
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The dependence of that parametrization reflects the following observations: i) it is constructed 
to describe the F2 data from BCDMS, NMC and HI which together cover a range of 4 orders 
of magnitude in x and Q"^; ii) for large x, F2is known to vanish like (1 — a;)^ with the parameter 
g close to 3 as determined by quark counting rules; iii) due to momentum conservation the 
integral over F2 is nearly independent of Q"^. Thus the rise of F2 with at low x must 
be compensated by a decrease of F2 with Q'^ at large x. As a consequence, at a; = Xg, F2 
is independent of Q"^. This demands that a polynomial term be included in front of the 
Q"^ dependent part which has been chosen to be (1 + e^yx). Since Xg in the flxed target 
experiments has been determined to be around 0.12, the parameter e should be close to —3; 

iv) the Q"^ dependence of F2 is expected to be logarithmic. The attempt to describe the Q"^ 
evolution over almost four orders of magnitude, from > 4 GeV^ to Q'^ < 2000 GeV^, 
requires a quadratic term in InQ"^. Fits were performed with and without that term present; 

v) flnally the introduction of a term in eq.7 is required which is independent of Q'^. 



a 


b 


c 


d 


e 
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g 


3.07 


0.75 


0.14 


-0.19 


-2.93 


-0.05 


3.65 



Table 5: Parameters of a phenomenological fit to the proton structure function data from 
tliis experiment combined witli _Ff from the NMC and the BCDMS experiments. The 
parametrization is valid for 4 GeV^ < < 2000 GeV^, lO'^ < x < 1 and < x -10^ GeV\ 

The result of the flt to the HI, NMC and BCDMS data is shown as functions of x and 
Q"^ in flgs. 6 and 7. The parameter values are quoted in table 5. In the flt the statistical 
and systematic errors of all quoted F2 values were added in quadrature and the relative 
normalizations were not allowed to vary. The flt provides a valid description of all data 
from the experiments considered here with a x^/rfo/ of 1.34. For the HI data alone the 
parametrization gives a x^/rfo/ of 0.85. The parameters g and e come out as expected to be 
close to ±3. From the result e = —2.93 one flnds that the slope of F2 with InQ"^ changes 
sign at Xg ~ 0.12. The obtained parameter / of eq.7 is small, i.e. the In"^ term amounts 
to about a 20% correction to the linear behaviour at Q"^ values between 20 and 50 GeV^. If 
the In"^ term in eq.7 is neglected the x^/dof increases from 1.34 to 1.55. An interesting 
result of the parametrization is the interplay between the flrst and the second term with an 
exponential x dependence. The second term dominates at low x. The value of the power d 
is found to be correlated with the presence of the In^ term. If the latter is neglected the 
flt suggests an x dependence at low x which is somewhat steeper, namely d = —0.23 instead 
ofd= -0.19. 

In flg.8 the F2 behaviour is illustrated for different Q"^ values as a function of W which is 
the invariant mass of the virtual photon-proton (j*p) system, 

W = ^Q^- {II X - 1) + ~ .Jcpfx (8) 

at low X. Here Mp is the proton mass. Since Q"^ /x = sy, the HERA experiments reach much 
larger W values than the previous DIS or real photoproduction experiments. F2 is displayed 
for the high statistics nominal vertex data and only for those Q'^ bins where x below 10~^ 
was reached, i.e. from 8.5 GeV^ to 35 GeV^. The rise of F2 at low x corresponds directly to 
a rise of F2 with W. This can be interpreted as a strong rise of the 7*p total cross section, 
since for a; << 1 and in the region where effects due to Z° boson exchange can be neglected, 

4 Tr^rv 

aUrp)^-7^F2{W,Q'). (9) 
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W / GeV 

Figure 8: Low Q"^ structure function data of HI plotted as functions of the invariant mass W 
of the 7*p system. The straight line is a ht F2(W, Q^) = 0.0058 W/GeV + 0.42 for x > 0.001 
which extends up to W c:^ 130 GeV where scales for the Q"^ range considered. Data for 
Q"^ > 80 GeV^ have a similar slope versus W but are below the lower Q"^ data. 

Three observations are made: i) the rise of with W is stronger than the one observed 
for W between 20 and 200 GeV of the total photoproduction cross section [23, 51], in which 
the photon is real. This behaviour for the total cross section of off-shell particle scattering 
was expected [52] and discussed for previous DIS experiments in [53]; ii) for W values below 
~ 130 GeV, corresponding to a; > 10~^, all the measured points cluster in a narrow band 
which is well reproduced by a straight line fit, i.e. F^iW, Q"^) scales in the Q"^ range considered. 
This observation is consistent with the recent ZEUS measurement [8]; iii) the extrapolation of 
the straight line fit (dashed line in fig. 8) into the higher W region, W > 130 GeV, reproduces 
the data at Q"^ > 20 GeV^ but the F2 data at lower Q"^ and x < 10~^ appear to deviate 
systematically from this linear behaviour. This trend is confirmed, though with less precision, 
by the lowest Q'^ measurements obtained with the shifted vertex data. 

The behaviour of at very low x requires still more precision data and extended coverage 
than were available to this analysis. A more theoretical discussion of the F2 data presented 
in this paper is deferred to a forthcoming analysis in the framework of perturbative Quantum 
Chromodynamics . 
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Figure 9: E'^, 9g and distributions for DIS events witli i]max < 1-8 and tlie-q^ax distribution 
for all DIS selected events. The Monte Carlo simulations are described in the text. 



7 DifFractive Contribution to F2(x,Q'^) 

In the previous HI measurement of F^ix^Q"^) [2], it has been observed that for about 6% of 
the deep inelastic events there was no significant energy deposition in the forward region. The 
angular acceptance of the LAr calorimeter is limited to polar angles ^ > 3° corresponding 
to a maximum measurable pseudo-rapidity, j] = — ln(tan|), of 3.6. The t] value of the 
most forward cluster [39] with energy greater than 400 MeV is defined to be i]max- The 
events with i]max < 1-8 amount to 6.3% of the low Q'^ data. A first HI analysis of the 
rapidity gap events collected in 1993 has been reported in [10]. There it was shown that 
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Figure 10: Measurement of the difFractive contribution F^^x^Q"^) to F^ix^Q"^) for Xpfp < 
10~^ and different Q"^ values witli the E method (closed circles) and the S method (open 
circles). The inner error bar is the statistical error. The full error bar represents the sta- 
tistical and systematic errors added in quadrature, not including the systematic error of the 
luminosity measurement of 4.5%. 



these events may be described by difFractive scattering processes where the virtual photon 
probes a colourless component of the proton, the Pomeron. The kinematics of deep inelastic 
difFractive scattering leads to a relationship [54] between the minimum size of the rapidity 
gap and Xpfp, the momentum fraction of the proton carried by the Pomeron, such that 

1 / \ 

A?7 = Vproton - Vmax > In with Xp/p = X ■ [1 + — ^ (10) 

where Mx is the invariant mass of the 7*iP system. This relation implies that a selection 
requiring i]max < 1-8 introduces a limit of Xpfp < 0.01 on the data [55]. The difFractive event 
sample contains both interactions with a scattered proton, and with a proton dissociative 
system. The analysis uses the same data sample and selection criteria as for the F^ix^Q"^) 
measurement with the additional event selection oi ri^ax < 1-8. 

Following [56] the difFractive contribution F^^x^Q"^) to F^ix^Q"^) for Xpjp < 0.01 can be 
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defined as 

D 2 1 

with the difFractive cross section and the kinematic factor k defined as in eq.6, setting 
i? = 0. Here t is the momentum transfer squared between the incident and scattered proton 
or proton dissociative system. To correct for the acceptance and for effects due to finite 
resolution, the RAPGAP Monte Carlo program [46] was used assuming a hard structure 
function for the Pomeron oc Xif]p(l — Xijp) with the same amount of quark and gluon induced 
events. Here Xijp is the fraction of the Pomeron momentum carried by the parton i which 
is assumed to interact with the virtual photon. A contribution from the elastic production 
of light vector mesons (/>(770), a;(783), (?!)(1020)), consistent with that observed in [10], was 
added to the RAPGAP Monte Carlo events [57] amounting to 10% of the selected rapidity 
gap sample. This simulation describes the data well as shown in fig. 9. Here the El,9e,y-E 
distributions for DIS events with i]max < 1-8 and the i]max distribution for all DIS events are 
shown for both data and Monte Carlo simulation. 

The contributions derived from the electron method and from the S method are 
shown in fig. 10. Both measurements agree well within the errors in the regions of overlap. 
As for the measurement of F2 the S method was taken for y < 0.15 and the electron method 
for larger y. The resulting values for are given in table 6. The systematic errors were 
calculated as for the inclusive F^ix^Q"^) measurement taking into account the uncertainties 
specific to the rapidity gap event analysis [55]. 
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Table 6: DifFractive contribution F^(x,Q''') to F^ix^Q"^) for Xpfp < 0.01. All points have an 
additional scale error of 4.5% due to the uncertainty of the luminosity measurement. 

By construction F^ decreases to zero at a; = 10~^. Comparison of the values of _F^(a;,Q^) 
with F^ix^Q"^) for x < 10~^ shows that the former contributes about 10% to the proton 
structure function. Therefore diffraction with Xpjj, < 10~^ cannot account for the steep 
increase of F2 with decreasing x. The Q'^ dependence of F^(x,Q''') at fixed x is fiat and 
not significantly different from the Q"^ dependence of the F^ix^Q"^). An analysis of diffractive 
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events in terms of structure functions will be given in more detail in a forthcoming publication 
of HI [58]. 

8 Summary 

A measurement has been presented of the proton structure function F2(x, Q"^) in deep inelastic 
electron-proton scattering at HERA with data taken in the running period of 1993. The 
integrated luminosity is 0.271 pb~^ which represents a tenfold increase in statistics compared 
to the first publication of HI. The structure function measurement includes data from 
different detector components and running configurations. Low Q"^ values are reached using 
data with the ep interaction vertex shifted in z from the nominal position. The data cover a 
kinematic range for Q"^ between 4.5 and 1600 GeV^ and x between 1.8 • 10"'* and 0.13. 

The F2 values presented are obtained using different methods to reconstruct the inclusive 
scattering kinematics. At high values of the scaling variable y > 0.15, due to its superior 
resolution, the electron method is used which is based on the scattered electron energy and 
angle. Lower y values are covered with the S method which combines electron and hadronic 
information to reduce radiative corrections and calibration errors. For comparison and cross- 
checks, F2 data obtained with the double angle method are also presented. 

The measured structure function has statistical and systematic errors considerably smaller 
than the previous result. For the first time, the HERA structure function measurement 
extends to the kinematical region of the high precision, fixed target experiments revealing a 
smooth transition between the BCDMS, NMC and HI results. 

A distinct rise is observed of the structure function with decreasing x at fixed Q'^. Around 
X ~ 10~^ the decrease of x by an order of magnitude amounts to a rise of F2 of about a factor 
two. This rise cannot be explained by the rapidity gap events because they contribute only 
about 10% of F2. They are used to measure for the first time the diffractive contribution F^ 
to the proton structure function for Xpfp smaller than 10~^. This contribution exhibits no 
significant Q'^ dependence. 

The observed Q'^ behaviour is consistent with the expected scaling violations, i.e. a weak 
rise of F2 with increasing Q"^ for x < 0.1. A parametrization is given of the proton structure 
function data from this experiment combined with the data from the NMC and BCDMS 
experiment describing F2(x,Q'^) over almost four orders of magnitude in x and Q"^. If F2 at 
low Q'^ is analyzed as a function of the virtual photon-proton mass W, deviations from a 
linear behaviour become visible for W > 130 GeV which deserve a more precise analysis with 
forthcoming data. 
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